Quantifying Quantum Correlations in Fermionic Systems using Witness Operators 
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We present a method to quantify quantum correlations in arbitrary systems of indistinguishable 
fermions using witness operators. The method associates the problem of finding the optimal entan- 
glement witness of a state with a class of problems known as semidefinite programs (SDPs), which 
can be solved efficiently with arbitrary accuracy. Based on these optimal witnesses, we introduce a 
measure of quantum correlations which has an interpretation analogous to the Generalized Robust- 
ness of entanglement. We also extend the notion of quantum discord to the case of indistinguishable 
fermions, and propose a geometric quantifier, which is compared to our entanglement measure. Our 
numerical results show a remarkable equivalence between the proposed Generalized Robustness and 
the Schliemann concurrence, which are equal for pure states. For mixed states, the Schliemann con- 
currence presents itself as an upper bound for the Generalized Robustness. The quantum discord is 
also found to be an upper bound for the entanglement. 

PACS numbers: 03.67.-a 
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I. INTRODUCTION 



The notion of entanglement, first noted by Einstein, 
Podolsky and Rosen is considered one of the main 
features of quantum mechanics, being subject of study 
in several areas recently 043- Thus, it is of paramount 
importance both the understanding and the quantifica- 
tion of entanglement in composite quantum systems, be- 
ing these one of the main challenges of modern quantum 
theory. 

Despite being a subject widely studied in systems of 
non-identical particles, or whose particles, though identi- 
cal, are well separated from each other, being thus distin- 
guishable, less attention was given to the case where this 
separation is very small, such that the overlap of their 
wave functions is no longer negligible. In this case we 
have to take into account the indistinguishability of the 
particles, being the space of quantum states restricted 
to symmetric or antisymmetric subspaces of the Hilbert 
space, depending on the bosonic or fermionic nature of 
the system. The study of entanglement for systems of 
indistinguishable particles has been, however, a subject 
of great controversy, leading to different approaches in its 
treatment. Among the most mentioned we have: quan- 
tum correlations [a], entanglement of modes |(|, entan- 
glement of particles Q. The concept of quantum corre- 
lations is based on the notion that the correlations gen- 
erated by mere (anti-)symmetrization of the state due 
to indistinguishability of their particles do not constitute 
truly as entanglement. We will analyse, in this paper, 
the quantum correlations in systems of indistinguishable 
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fermions. 

The entanglement in two-fermion sates with a four- 
dimensional single-particle space (ACH^ ® 'H^ 4 ]Y) can 
be characterized by the Schliemann concurrence [5[ . For 
pure states in arbitrary dimension, one can use the von 
Neumann entropy of the reduced single-particle density 
matrix as a measure of entanglement @. However, the 
problem of quantification or even detection of entangle- 
ment in the general case is still open A useful 
concept is that of entanglement witness [12[ : a Hermitian 
operator with non-negative expectation value for all sep- 
arable states, but which can have a negative expectation 
value for an entangled state. We will focus on optimal 
entanglement witnesses (OEW), that can be used not 
only to witness the entanglement, but also to quantify it 
jFfl - ll5l ]. In this paper, we will see how to determine such 
OEWs, and especially how to use them to quantify the 
quantum correlations in fermionic systems. We will also 
confront our measure of entanglement with the quantum 
discord EMI. In order to do that, we will define a 
quantum discord for fermionic particles. Previous stud- 
ies of quantum discord in systems of indistinguishable 
particles, like in [2(| , employ the notion of correlation of 
modes [6(, which are distinguishable. 

This paper is organized as follows. In Sec|TT]we recall 
some known concepts and tools of the theory of entan- 
glement, as the witnessed entanglement and the quan- 
tum correlations, which will be essential to the devel- 
opment of the ideas throughout the article. In Sub- 
scc lll A[ we recall the definition of optimal entanglement 
witness, and briefly discuss its use as an entanglement 
quantifier. In Sub-sec fll Bl we recall the ideas associated 
with the concept of quantum correlations, as well as the 
Schliemann concurrence. In Sub-sec llI CI we revise the 
concept of quantum discord for distinguishable particles. 
In Sec lIIII we introduce our method for quantifying the 
quantum correlations in fermionic systems using witness 
operators, and define the Fermionic Generalized Robust- 
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ness. We also extend the notion of quantum discord 
to fermions, such that it takes into account the parti- 
cles' indistinguishability, and introduce the Fermionic 
Geometric Discord. In Scc lIVI we show numerical re- 
sults, comparing the Fermionic Generalized Robustness, 
the Schliemann concurrence, and the Fermionic Geomet- 
ric Discord. We finish the illustration of the method 
with a beautiful quantum phase diagram yielded by the 
five-partite Fermionic Generalized Robustness in the Ex- 
tended Hubbard Model. We conclude in Sec.V. 



II. PRELIMINARY CONCEPTS 



EqU using M = {W G W | W < Z}Q] 
identity operator; in other words, 



where I is the 



Tt g {p) = max{0, 



min Tr(Wp)). (4) 

{wew\w<i} 



The construction of entanglement witnesses is a hard 
problem. In an interesting method proposed by Brandao 
and Vianna hlj, the optimization of entanglement wit- 
nesses is cast as a robust semidefinite program (RSDP). 
Despite RSDP is computationally intractable, it is pos- 
sible to perform a probabilistic relaxation turning it into 
a semidefinite program(SDP), which can be solved effi- 
ciently m. 



We will see, in this section, some familiar concepts from 
the theory of entanglement, namely, witnessed entangle- 
ment fl3j-|l5| , quantum correlations Q and quantum dis- 
cord [16l4l9j . A reader already familiarized with such 
concepts might skip to the next section of the article. 



A. Witnessed Entanglement 

Entanglement witnesses are Hermitian operators (ob- 
servables - W) whose expectation values contain infor- 
mation about the entanglement of quantum states. The 
operator W is an entanglement witness for a given en- 
tangled quantum state p if the following conditions are 
satisfied [12( : its expectation value is negative for the par- 
ticular entangled quantum state {Tr{Wp) < 0), while 
it is non-negative on the set of separable states (S) 
(V<7 G S, Tr(Wa) > 0). We are particularly interested 
in optimal entanglement witnesses. W op t is the OEW for 
the state p if 



Tr(W optP ) 



min Tr(Wp), 
weM 



(1) 



where A4 represents a compact subset of the set of en- 
tanglement witnesses W. 

OEWs can be used to quantify entanglement. Such 
quantification is related to the choice of the set M. , where 
different sets will determine different quantifiers [13( ■ We 
can define these quantifiers by: 



E{p) = max(0, 



min Tr(Wp)). 
WeM 



(2) 



An example of a quantifier that can be calculated u sing 
OEW is the Generalized Robustness of entanglement [22[ 
( IZg(p)), which is defined as the minimum required mix- 
ture such that a separable state is obtained. Precisely, it 
is the minimum value of s such that 
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(3) 



be a separable state, where tp can be any state. We know 
that the Generalized Robustness can be calculated from 



B. Quantum Correlations 

The space state for indistinguishable fermions is anti- 
symmetric under permutation of particles. In this case, 
it is convenient to use the second quantization formalism, 
in order to deal with the antisymmetric states in the Fock 
space. Accordingly we introduce an algebra of operators 
which satisfy the following anti-commutation relations: 



{//,/]} = {/i,/;} = 0, {/i,/j} = % 



(5) 



where /• and fi are the fermionic creation and annihila- 
tion operators, respectively, so that their application on 
the vacuum state (|0)) creates/annihilates a fcrmion in 
state "i". The vacuum state is defined so that fi |0) = 0. 

An immediate and disturbing consequence of the anti- 
symmetric structure of the state space can be seen even 
in the simplest example of a two-fermion system, which, 
if analysed in the usual way, will always be considered 
entangled. We must therefore rethink the way entangle- 
ment is calculated for systems of indistinguishable parti- 
cles, as well as its physical interpretation. 

In the case which the identical particles are localized 
in distinct laboratories and independently prepared, it is 
natural to think that the entanglement calculated in the 
usual way should not have any relevant physical meaning; 
or rather, "no quantum prediction, referring to an atom 
located in our laboratory, is affected by the mere presence 
of similar atoms in remote parts of universe" [24j . 

We are interested in the case of identical particles that 
are sufficiently close together such that the overlap be- 
tween their wave functions is no longer negligible, and 
therefore they are indistinguishable. Fermionic systems 
of this kind can be described using Slater determinants 
3. Consider, for example, a two-fermion state repre- 
sented by a single Slater determinant, namely, 

H') = 4=(M ® \x) - \x) ® m = flfl |0) , (6) 

where \4>) and |%) correspond to orthonormal wave func- 
tions (spin-orbitals). It is easy to sec, in this simple 
case, that the anti-symmetrization of coordinates intro- 
duces correlations between the fermions, namely, the well 
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known exchange contributions from the Hartrcc-Fock the- 
ory. On the other hand, a single Slater determinant is so- 
lution of a one-particle Schrddingcr equation and, there- 
fore, can have no quantum correlation between the par- 
ticles [ij . Considering states described by more than one 
Slater determinant introduces additional correlations be- 
yond the exchange contribution. We will then interpret 
such additional correlations as the analog of quantum en- 
tanglement in systems of distinguishable particles, calling 
them as fermionic entanglement || . 

A measure of fermionic entanglement was proposed 
in Q as the analogous of Wootters concurrence f25j . 
Notwithstanding, such measure, called Schliemann con- 
currence (Cs), is valid only for two-fermion states with a 
four-dimensional single-particle Hilbert space {A(H^ ® 
W^)), i.e. the antisymmetric space of lowest dimension 
where can exist quantum correlated states. 

In order to define the Schliemann concurrence, we have 
to introduce some operators. Let U p h be the operator of 
particle- hole transformation: 

d 

Uphf\i4h = & u ? h i°) = n ft i°) > ( ? ) 

i=l 

being d the single-particle Hilbert space dimension. Sim- 
ilarly, define K. as the anti-linear operator of complex- 
conjugation, satisfying the following relations: 

KflK = fl KfiK = fi, K|0> = |0>. (8) 

Thus, given the operator T> = JCU p h, called operator of 
dualisation, and the dual states p = PpD -1 , we have that 
the Schliemann concurrence for states p e .4("H (4) ®H (4) ) 
is given by 

C's(p) = max(0, A 6 - A 5 - A 4 - A 3 - A 2 - Ai), (9) 

where X^s arc, in descending order of magnitude, the 
square roots of the singular values of the matrix R = pp. 

C. Quantum Discord 

The total amount of correlations (classical and quan- 
tum) carried in a bipartite system is quantified by the 
mutual information: 

I(A:B) = S(A) + S(B)-S(A,B), (10) 

where S(X) is the Shannon (von Neumann) entropy for 
the classical (quantum) system. In the classical case, the 
mutual information has another interpretation: it mea- 
sures the decrease of ignorance about the subsystem A, 
when B is known. In this case, an equivalent expression 
for the mutual information reads: 

J(A : B) = S(A) — S(A\B), (11) 

where 

S(A\B) = Y,PbS(A\b) (12) 

b 



is the conditional entropy on A given B. In quantum 
systems this equivalence is not always true. The dis- 
agreement between these expressions is called quantum 
discord, and can be used to quantify quantum correla- 
tions [161 ] : 

D(A : B) = I (A : B) - J {A : B). (13) 

However, there is an ambiguity in J(A : B) due to the 
freedom of choice for the measurement operators, where 
for each set of measurements {n^js on B, the conditional 
entropy value may be different and consequently also the 
mutual information. The quantum discord is thus rede- 
fined by minimizing its value over all measurement oper- 
ators. The mutual information is thus properly defined 
by: 

J (A : B) = S{A) - max S(A\B), (14) 

and it is known to quantify classical correlations [l7j . 
Thus, the quantum discord measures the difference be- 
tween the total and classical correlations on the sys- 
tem, therefore being a quantifier of quantum correlations. 
Note however that quantum discord is not a measure 
of entanglement, and quantum separable states usually 
have non-zero discord. The quantum states with null 
discord are those which are a mere encoding of classical 
statistical distributions, and can be written in the form: 

£=5>^>< e *i®i£><£i- ( 15 ) 

ij 

where \ei) and \ fj) are orthonormal basis. 

Given the set (f2) of states with zero-discord, formed 
by the states above, it is clear that the minimal distance 
between a state p and this set can be used to quantify 
the quantum discord [l8j : 

where \A\ P is the Schatten p — norm. A usual measure 
of quantum discord is based on the 2-norm, or Hilbert- 
Schmidt norm, proposed by Dakic et al. [19( . 

III. QUANTIFYING QUANTUM 
CORRELATIONS USING WITNESS 
OPERATORS 

In this section we will present our method for quan- 
tifying fermionic entanglement using witness operators. 
After defining the set of separable states pertinent to our 
case, i.e. the fermionic states without entanglement, we 
will introduce a method for determining OEWs. In par- 
ticular, we will see which constraints to impose on the 
set A4 of witnesses (Eq[2]) in order to obtain a quanti- 
fier analogous to the Generalized Robustness (Eqf?]). We 
also will define a notion of quantum discord for fermions 
and confront it to entanglement. 
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We know that states without entanglement are those 
that can be described by a single Slater determinant, 
or a convex mixture of them. Consider as the 
Fock space of n indistinguishable fermions sharing a 
(i-dimensional single-particle Hilbert space. We have 
then the following definition of "separable" states: 

State with no fermionic entanglement (separa- 
ble): A state a £ B(J-^) has no fermionic entanglement 
if it can be decomposed as 



E 



it 

Pi a x 



where a 



T 



d 



• •<ho)(o| £> = i, 

(17) 

// , and {/, } is an ortho normal basis 



of fermionic creation operators for the space of a single 
fermion (F?). Of course, the states defined by Eq flTl are 
not separable in the usual mathematical sense, mean- 
ing that they are product states or convex mixtures of 
it. But we will insist in referring to them as separable, 
for they are just anti-symmetrization of the usual distin- 
guishable separable states. Entanglement, in the case of 
distinguishable particles, is defined in opposition to sep- 
arability, i.e., an entangled state is that one which is not 
separable. We want to keep this notion. 

It is interesting to note that, as in the case of distin- 
guishable particles, the set of separable states is invariant 
under local operations, taking now into account that the 
local operations must be symmetric, due to the indistin- 
guishability of the particles. Let $ be a local symmetric 
operation (LSO), i.e., an operation that respects the Pauli 
exclusion principle and does not involve any interaction 
between particles. An LSO can be written as: 



Hp) = £(Mi 



■®M i )p(M}®M i 



rt, 



■®M}) 
(18) 

where M, is a linear operator acting on the Hilbert space 
of a single particle. Given a fermionic separable pure 
state (i.e. a single Slater determinant) \tp sep ) = v4(|ei) ® 
•••<£> |e n ))j where A is the anti-symmetrization operator, 
{|ej)} is an orthonormal basis, and noting that [$,-4] = 
0, we see that 



(Mi ® • • • ® Mi) IV'Sep) 



= (Mi® •••(8 


Mi)A\e x ---e n ) 


= A (Mi ® ■ ■ • 


<g> Mi) |ei • • • e„) 


= A |(Mei)- 


• • {Mien)) 


= A\e' 1 ---e' n ) 


(19) 


= Wsep) 





and such result clearly extends to mixed states. Sum- 
marizing, given a separable state a £ S, we have that 
^lso( <t ) € <5, indicating that in order to have quantum 
entanglement, the particles must interact by means of 
some global operation. 



Now we adapt Brandao and Vianna's [141 ] technique in 
order to obtain a new algorithm to determine OEWs for 



indistinguishable fermions in the Fock space. The new 
method can be enunciated as follows. 

Determination of OEW using RSDP: A 
fermionic state p £ B(T^) is entangled if and only if 
the optimal value of the following RSDP is negative: 

minimize Tr(W p) subject to 



d d 



E ■■■ E E ••• E (C-i -«4'i-C x 

l,i_l=l ii = lj'i = l j„_ 1 = l 

Win-l"l-lil-J»-l)>0, 

Vc- £ C, l < k < (n - l), l < i < d, 
AWA^ = W, 
W <A, 

(20) 



where d is the dimension of the single particle Hilbert 
space, {//} is an orthonormal basis of fermionic creation 
operators, A is the anti-symmetrization operator, and 
W^-iiii-in-i = fi n - l ---fi l Wf] l ---fl_ l £ B(F?) 
is an operator acting on the space of one fermion. The 
notation W < A means that (A — W) > is a positive 
semidefinite operator. If p is entangled, the operator W 
that minimizes the problem corresponds to the OEW of 
P- 

Proof: It is known that a state is entangled if and only 
if there exists a witness operator W such that Tr(Wp) < 
and Tr(Wa) > for every separable state a. Consider 
a general separable state as given by EqlTT] The semi- 
positivity condition Tr(Wo~) > is equivalent to: 



| a„a„_i ■ ■ -ai W a\ ■ ■ -a^jO* |0) > 0, 



(21) 



for all a\ £ B(Tf). Note however that to satisfy Eq|2"T1 
it is sufficient that the operator a„_i ■ ■ • a± W a\ ■ ■ ■ a\ l _ 1 
be positive semidefinite. Thus follows directly that the 
operator W satisfying the problem in Eq |2"01 corresponds 
to an optimal entanglement witness. 



The RSDP given above is solved by means of prob 
abilistic relaxations it terms of SDPs, as done in 
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where the set of infinite constraints is exchanged by a 
finite sample. Thus the witness operator obtained is 
such that satisfy most of the constraints in Eql^Ol The 
small probability (e) that a constraint be violated (i.e. 
Tr(Wa) < 0) diminishes as the size of the sample of 
constraints increases. 

The constraint AWA^ — W restricts the operator 
to the space of antisymmetric entanglement witnesses 
(W(J^) = AWA'). The other constraint, W < A, 
follows directly from the anti-symmetrization of Eq|4l 
and implies that the OEW corresponds to the anti- 
symmetrized version of the Generalized Robustness, 
namely, 

TZ^(p) =max(0,- min Tr(Wp)). (22) 

9 M = {W£W(J r £) | W<A} 
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FIG. 2: Fermionic Generalized Robustness IZg for the families 
of two-fermion states defined in Eg 1261 on the space A{w- 2L ^ ® 



distributions. Although these states do not have any 
kind of quantum correlations, they cannot be treated 
like classical probability distributions, since they respect 
quantum rules: like the Pauli exclusion principle. 

Our proposed measure for the quantum discord in 
fermionic states will be a geometric measure like Eq ll6l 
Given a fermionic state p € A(7i^<S>'H^), the Fermionic 
Geometric Discord is given by, 



FIG. 



1: Fermionic Generalized Robustness lZ g versus Schlie- 



mann Concurrence Cs, for random fermionic states uniformly 
distributed according to the Haar measure. (TOP) The two 
entanglement measures are equal for pure states. The small 
dispersion seen in the top panel is due to numerical impreci- 
sion in the calculation of IZf in the region of very low entan- 
glement. (BOTTOM) In the case of mixed states, the Cs is 
an upper bound to IZg ■ The continuous line in the bottom 
panel corresponds to the straight line Cs = TZf , and is just a 
guide to the eye. 



TZg(p) measures the minimum required mixture with a 
fermionic state such that all the entanglement of p is 
washed out. In other words, the Generalized Robustness 
is the minimum value of s such that 



3<p f 



1 



(23) 



be 



any 



be a separable state (EqfT7|), where <pj 
fermionic state. 

Fermionic Geometric Discord: Abiding by the 
notion that mere anti-symmetrization does not gener- 
ate any kind of quantum correlation, we are led to the 
following definition of fermionic states without quantum 
discord: 



(24) 



where A is the anti-symmetrization operator, and £ are 
the states in EqfTSl which encode classical probability 



v f(p) = pip \€-p\u 



(25) 



where fi^ = AHA' is the set of zero-discord anti- 
symmetric states (Eq[ 



IV. NUMERICAL RESULTS 

In this section we will illustrate our method. We 
start by investigating bipartite entanglement in the space 
A(H^ ®H.^), which has the smallest dimension allow- 
ing for quantum correlations in fermionic systems. In 
this case, we can compare the Fermionic Generalized Ro- 
bustness {Jig - Eq |22|) with the Schliemann concurrence 
(Cs - Eq|9|). Then we investigate the bipartite entan- 
glement in a one-parameter family of states in the space 
-4(H (2i) <g)'H (:2i) ), with L going from 2 to 4. We also com- 
pare our Fermionic Geometric Discord (Df - Eq |25l) with 
the Fermionic Generalized Robustness for another one- 
parameter family of states in the space A^H^ (£> H^). 
We finish with calculations of multipartite entanglement 
in the Extended Hubbard Model (EHM), where 7l£ can 
characterize quantum phase transitions. 

In Fig.l, we plot the Fermionic Generalized Robustness 
against the Schliemann concurrence for states generated 
randomly in the space AiH^ ®'H^). It is remarkable 
that Cs and IZ^ are one and the same for pure states, 
and Cs is an upper bound for TZ^ in the case of mixed 
states. Recall that, in the case of two distinguishable 
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FIG. 3: Formionic Generalized Robustness IZg , and 
Fermionic Geometric Discord 2?/ for the families of two- 
fermion states denned in Eq |26l (top) and Eg 1271 (bottom), 
on the space A(H {4) ® H (i) ). 



qubits, the Wootters Concurrence [25( and the General- 
ized Robustness [22| also keep this same relation. 

Now we consider the following one-parameter family of 
states, in the space A(n {2L) ® H^): 



P = f(0>+f(l/2)p max +f(l)A, f(x)=Ae A 2 , 

(26) 

where a = fjfj |0) (i ^ j) is a pure state with just a sin- 
gle Slater determinant; p m ax is the maximally entangled 
pure state of singlet type, i.e. the one with spin quantum 
numbers S = S z = 0; and A is the anti-symmetrizer, 
which corresponds to the identity operator in the an- 
tisymmetric space, p controls the entanglement of the 
state, A is chosen such that the state is normalized, and 
A = 0.1826. In Fig. 2, the Fermionic Generalized Ro- 
bustness is plotted against p, for spaces of different di- 
mensions. We see that behaves correctly, showing 
that the mixed state has much entanglement when the 
contribution of the singlet is large, and it has low, or 
none, entanglement when the contribution of cither a or 
A is large. 

Consider another one-parameter family of states, in 
the space A(H {4) ® H {i) ): 



P = (1 -p)A + pp v 



< p < 1, 



(27) 



In Fig. 3, we confront the Fermionic Generalized Robust- 
ness and the Fermionic Geometric Discord, for the two 
families of states (Eq.26, Eq.27) defined above. In the 



FIG. 4: Five-partite Fermionic Generalized Robustness in the 
Extended Hubbard Model (Eq.19). The multipartite entan- 
glement works as an order parameter and characterizes quan- 
tum phase transitions, in this case, with three distinct phases. 



two cases, we see that the discord is always an upper 
bound for the entanglement, and states without entan- 
glement can have a non-null discord, which is particu- 
larly dramatic in the second family. In the first family of 
states, the functional forms for the discord and entangle- 
ment are very similar, whereas they are very different in 
the second family. Note, in the second family, the abrupt 
vanishing of entanglement (a discontinuous derivative for 
p = 0.8), while the discord shows an asymptotic behav- 
ior. 

To conclude, we illustrate the calculation of multipar- 
tite entanglement of fermions interacting according to the 
Extended Hubbard Model [2l| , denned by the Hamilto- 
nian: 

H E HM = ~th J2j = l, < j= J (,i(ftafj+l,'r + /j+l, CT /.J>) + 

u Ej=i nftUjx + V £\=i n 3 n j+1 , 

. (28) 

where U and V define the on-site and nearest-neighbor 
Coulomb interactions, th controls hopping between sites, 
L is the number of sites, and rij a is the particle number 
operator on site j with spin a. Fig. 4 shows the five- 
partite Fermionic Generalized Robustness of the ground 
state as a function of U /th and V/th, for the case of five 
fermions and five sites. The ground state is obtained by 
numerical diagonalization of the Hamiltonian. One can 
see that the Fermionic Generalized Robustness charac- 
terizes three distinct regions, corresponding exactly to 
the three distinct phases provided by the known phase 
diagram of the model, namely charge-density wave (up), 
spin-density wave (right) and phase separation (bottom) 
[2l|. It is interesting to note that five fermions in five 
sites, with periodic boundary conditions, is the small- 
est size of the system which presents such phase transi- 
tions. The occurrence of these phase transitions is depen- 
dent on long range interactions. Surprisingly, the second- 
neighbor interactions present in our five fermions model 
is already long range enough for the onset of phase tran- 
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sitions. We performed a calculation with four fermions 
in four sites, and the resulting figure is a uninteresting 
flat surface for the entanglement. 

V. CONCLUSION 

In summary, we presented a method to quantify quan- 
tum correlations in systems of fcrmionic indistinguishable 
particles. The method is based on the use of optimal en- 
tanglement witnesses, which can be calculated with arbi- 
trary precision by means of SDPs. In particular, we ob- 
tained the Generalized Robustness for fcrmionic systems 
(IZg ), and numerically showed its relation to the Schlie- 
mann concurrence. We also introduced the Fermionic 
Geometric Discord (2?/), and observed that it is an up- 



per bound for the fermionic entanglement. However, the 
physical meaning of quantum discord for fermionic sys- 
tems needs to be better understood. It is comforting to 
know that the quantum discord for a single Slater de- 
terminant, or for the fermionic maximally mixed state, 
is null, but the nature of quantum correlations conveyed 
by a convex mixture of Slater determinants is still ob- 
scure to us, and will be investigated in future works. 
Finally, we used the five-partite Fermionic Generalized 
Robustness to characterize quantum phase transitions in 
the Extended Hubbard Model, showing the utility of en- 
tanglement as a quantum order parameter. 
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